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Abstract 



I. INTRODUCTION 



The aim of this paper is to show that in the special case of stationary asymptotically 
flat spacetimes there are invariantly defined expressions for both the energy density of the 
gravitational field and for the angular momentum density. Furthermore, these expressions 



However the gravitational expressions do involve the gradiant of a conformal factor which is 
defined non-locally using a geometrical structure depending on the Killing vector. We hope 
that our expressions can be generalized to slowly varying systems in general relativity. We do 
not aim to introduce any new generally defined concept of energy or quasi-local energy such 
as those of Penrose, Hawking or Hayward [1]. We are concerned to isolate the gravitational 
potential energy in the stationary case and reexpress it in terms of 'fields' whereas those 
quasi-local expressions contain material energy as well as gravitational energy and go over 
to the ADM mass at infinity. 

Misner, Thorne and Wheeler (1973) [2], hereafter MTW, deny the existence of localised 
gravitational field energy-density (with somewhat strident rhetoric). Nevertheless they give 
an expression for the localised gravitational energy in the exceptional case of spherical sym- 
metry (quoting Misner & Sharp (1964) [3] to include time dependence). Katz (2005) [4] has 
recently given an expression for it in the more general case of conformastationary metrics. 
As our expression (which agrees with his) differs from that of MTW even for the spherically 
symmetric case, we shall start by considering their special case in this introduction. 

Landau & Lifshitz (1966) [5] show that the general spherically symmetric metric can be 
put in the form (setting c 2 = 1 and writing £ 2 for their e u ) 



are the gravitational analogues of (E 2 +B 2 



)/(8n) and [fx (ExB)]/ (Att) for electromagnetism. 




where m(r, t) tends to = GM^ the total mass at infinity times G, and 



r = (sin 9 cos tp, sin 9 sin ip, cos 9) , (dr) 2 = d9 2 + sin 2 9dtp 2 . 



The relevant component of Einstein's equations reads 




so 




(1) 



2 



Tq is the energy density of the matter in the rest frame (r, 9, 0) constant. When r = 2m(r, t) 
at some radius r*(i) a central black hole of that radius and mass M*(t) = r*/2G = m*/G 
is present. Such a case is considered in [6] but here we shall confine ourselves to global 
complete spacetimes. 

To see exactly what equation ([I]) means consider for example a fluid of dust. Then 
T^ v = pu^u u where p is the density of the fluid in its rest frame and is its four velocity 
which is of the form (u°, u r , 0, 0) by symmetry. 

Let w M be the 4-velocity of static observers with (r, 9, 0) constants, so these are unit vec- 
tors and in their coordinates have only one component = 0, 0, 0) , w^ = (£, 0, 0, 0). 
Since the fluid moves radially in our frame, they see it to have a greater density due to 
Lorentz contraction. This increased density of rest mass is pu^w^ = pu.w and the mechan- 
ical energy density in any frame now including kinetic energy is p ( u - w ) . This is the time 

y l—v 2 /c 2 

component of 

pu.wu" = T;< (2) 

which is the flow vector of mechanical energy and momentum seen by the static observers. 
The mechanical energy is its flux through a space-like hypersurface such as t = constant. 
The part of that flux through t = constant that lies within r = constant will be 

E M (r) = T°w»V=9 d 3 x = £ T °4vrr 2 - ^M^) ^ dr . ( 3 ) 

Em{t) is a special case of the energy seen by observers moving orthogonally to a hypersurface 
like that geometrically defined by Wald [12]. Although for illustrative clarity we considered 
only dust above, this is still the mechanical energy within r whatever constitutes the T^ v . 
Indeed for a perfect fluid it includes the internal energy as well as rest mass and kinetic 
energy, 

_ p + pv 2 

where v is the 3-velocity of the fluid in their frame. 

A comparison between equations (CQ) and ([3]) illustrates that (CQ) is seductively like the 
classical relationship between density and mass but in fact conceals all the complications 
beneath a cloak of apparent clarity. 

Firstly Tq is not the density in these coordinates unless v is zero; secondly 47rr 2 dr is not 
the volume element which is rather 

47rr 2 [l - 2m(r, t)/r]- 1/2 dr = dV . (5) 



Outside the matter m becomes constant and gives G times the total mass-energy seen 
from infinity. MTW, quoting Misner & Sharp [2], argue that M(r,t) = G~ 1 m(r, t) is the 
total mass-energy within r at time t and is only less than EM{r,t) because of Eo(r,t) the 
gravitational energy which is negative. Indeed, since Em (equation (j3J)) contains the rest- 
mass energy, the kinetic energy and the internal energy within the matter, the difference 



E M (r,t) -G- L m(r,t) 



1 



2m(r, t) 



-1/2 



1 [ T °47rr 2 cfr 



-E G (r,t) (6) 



must be due to gravitation. 

This would lead to a gravitational energy density of 



dE G (r) 
dV 



1 



1 



2m(r, t) \ 



1 o 



1 



1 



2m(r, t) \ 



AirGr 7 



dm 
dr 



which is only non-zero inside matter. 

To see that this deduction might not be water-tight, we turn to the electrostatic analogue 
in flat space. The electrical energy of a spherical charge distribution, Q(r) being within r, 
can be calculated by starting at the centre and imagining the distribution to be built up by 
adding shells of charge consecutively. The shell with charge dQ is added when the potential 
at r is Q(r)/r so the electrical energy up to r is 

f r Q(r) dQ 
Jo 



E Pm (r) 



-dr. 



r dr 

As further shells are added we merely extend the upper bound of the integral so the part 
within r does not change. However it would be wrong to deduce that the energy density in 
the electric field at r in flat space is 

1 dE em (r) _ Q dQ 
47rr 2 dr 4irr 3 dr 

An alternative is found from the formula E* m {r) = J Q r ^(f)(dQ/dr)dr where is the electrical 
potential = J^°(Q/r 2 )dr. This formula takes account of the fact that the electrical 
potential changes from Q/r as outer shells are added. However it is also wrong to imagine 
that the electrical energy density is {Anr 2 )- l (dEl m / dr) = (^r 2 )- l (j)(dQ / dr). Notice that 
both of the above formulae only give contributions from within the charge distribution Q(r) 
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whereas the true answer due to Maxwell is (87r) _1 J E 2 dV = (87r) _1 J (Q/r 2 ) 2 dV . Of course 
all three expressions integrate over space to the same total electrical energy, what is in 
dispute is the distribution of that energy. This electrical analogy suggests that it might be 
appropriate to evaluate to infinity and integrate the result by parts. Then it might be 
possible to express the result as the volume integral of a perfect square which might be 
interpreted as in the electrical case. 

Doing this we write x = 2m(r,t)/r and use ([!]) 



-GE, 



G 



poo 
JO 



-1/2 



l]dm/dr 



If 00 f dx\ 

dr = - J [(1 - xY 1/2 - 1] ( x + r—J dr 



[(1 



x 



,-1/2 



d 



l]x + r— [2-2(1 
or 



x 



) 1/2 -x] \dr . 



Integrating by parts and remembering that x — > as r — ► oo the two terms recombine to 
give us just such a perfect square integrated over the volume dV given by (jSJ) viz 

2 



- E, 



G 



8nGr 2 



1 



1 - 



2m(r, t) \ 



dV 



8ttG 



F 2 dV 



(7) 



where we have introduced a gravitational field strength F given by 



2m\ 
r J 



1/2 



~9rr) 



-1/2 



9$ 

dr 



where $ is a gravitational potential. Equation ([7]) may be seen as the gravitational analogue 
of Maxwell's J E 2 dV/ (8tt) but in the special case of spherical symmetry. For the geometrical 
interpretation of $ see below. An expression equivalent to (J7J), but expressed in isotropic 
coordinates was given in [4]. Denoting by — > values at large r, 



F r = -d§/dr 



2m 
r 



-1/2 



m 



Evidently 



oo ^ 

r 



1 



2m(r, t) 



-1/2 



1 > dr 



m 
r 



(8) 



What physical meaning should be ascribed to this gravitational field and what is the 
source of $? 



V 2 <T> ' f I 2m V /2 9 
r 2 \ r j dr 



2m\ 
r J 



1/2 



,5$ 

dr 



1 dm 1 



2m 

r 



r J 
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so 



- V •!• - I («T« - F 2 ) = \ [kT ° - (V$) 2 ] , 



2 v ~~ u " 7 2 L '~" u v ' _/ J ' ^ 
where k = 8nG. This shows that the negative gravitational field energy acts alongside 
the positive Tg as a source for $. In this last form it is natural to generalise (Q out of 
spherical symmetry and this is indeed what has already been done in [4] for conformastatic 
and conformastationary systems. 

We now turn to the interpretation of $. All spherical spaces have conformally flat spatial 
metrics. If we rewrite our spatial metric in isotropic form, 

da 2 = e 2A (dr 2 + r 2 dr 2 ) , 



then, at constant t, 



and 



so 



2m 
r 



r = e A r 



-1/2 



- 1 



dr 
r 



d In r — d In r = d In 



Hence, taking 



1 at infinity, we find, cf equation 



exp < — 



r / 



(10) 



So $ is A and e 2 * is —g\\ in isotropic coordinates. For time dependent spherical metrics 
r is geometrically defined as a real radius and the surfaces of constant t are geometrically 
defined as those symmetric cuts of spacetime orthogonal to r. 
If we make the conformal transformation da 2 = e~ 2 ^da 2 then 



da 2 = (dr 2 + r 2 dr 2 ). 



(11) 



Thus the interpretation of $ for the spherically symmetric case is that e -2 * is the con- 
formal factor of the transformation which makes the spatial metric flat. In general, spatial 
metrics of general spaces cannot be made flat, most of those that can have already been 
explored in [4]. It is worth noticing that the energy (|7j) or (J3j) depends solely on the spatial 
metric, £ is not involved. 



Although we can not generally make a conformal transformation to give a flat 3-space, 
we can in practice make a conformal transformation to a 3-space which has no scalar 3- 
curvature. This turns out to be a crucial step in determining the invariant energy density 
of general stationary gravitational fields. 

In Katz, Lynden-Bell & Bicak (2006) [6] this technique is employed to give energy 
densities according to both static and hypersurface orthogonal observers, such as ZAMOs 
(Bardeen (1970) [7]). 

Here we give a more physical exposition of that work and extend it to include angular 
momentum density. This leads us to a new and different expression for gravitational field 
energy V. 

For stationary space times we write our metric in the alternative forms 

ds 2 = £ 2 (dt — Akdx k ) 2 — r )kidx k dx l = g^dx^dx 
= ( 2 dt 2 - %i(dx k - W k dt)(dx l - W l dt) 

where Greek indices run from to 3 and latin ones from 1 to 3. Evidently 

Ak = -gok/goo and £ 2 A k = -gok = -IkiW 1 = -£ k (13) 
where £ M is the 'stationary' Killing vector. We define 

7 = det(7 fc i) 

and the equivalent for 7. 

Landau and Lifshitz [5] show that for a given spacetime there are many stationary metrics 
of the form (Tl2l) since an arbitrary function of the x k can be added to t without destroying the 
stationary form of the metric. Thus there are many different slicings of the given spacetime 
into time and space. However, as discussed in [6], we can get a unique geometrically defined 
slicing by demanding that it be maximal in the sense that the trace of the the external 
curvature, K, of the constant time slices be zero. This choice of time slicing is clearly a 
good one in that it picks out Boyer-Lindquist time in Kerr spacetime. We shall thereafter 
make this choice of the time coordinate so that the space on each constant time slice is a 
well defined geometrical concept. 
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II. GRAVITATIONAL FIELD ENERGY DENSITIES 



In [4] it was realised that the definition of mechanical energy used in spherical symmetry 
by MTW [1] could be extended to any stationary space-time. If w M = then the density 
of rest mass in a dust fluid seen by a static observer is -/j^p* = p( u - w ) and the energy 



density on any spacelike surface element dH^ is p(u.w)u fl ^—gdTj fl = T^Wy^J—gdU^ [12 ]. 
If we have a more complicated T^ u such as a gas or a plasma this last expression still 
measures the mechanical energy density as seen by the static observers but it now includes 
internal energy (including rest mass energy) and the kinetic energy of motion relative to the 
observers. If we sum these local contributions we get the total mechanical energy with no 
contribution from the gravitational binding energy 

E M = J T»w u V^gdJ:^ = J T ° dV, where dV = ^d 3 x. (14) 

We notice that this expression agrees with that given for spherical symmetry ([3]) and ([6]). 
However, this is not generally the same as Wald's expression because £ M and are not 
generally hypersurface orthogonal, so are not normal to the hypersurface. When black holes 
are absent we may now define the total gravitational energy of any stationary space-time 
by Eg = M — Em, where M is the total mass. We get an interesting expression for Eg by 
using Einstein's equations to transform Tq into 'field' quantities via integrations by parts 
analogously to Maxwell's treatment in electrodynamics. 

We write V for the vector operator in the 3-space with the 7^ metric and then 
follow Lynden-Bell & Nouri-Zonoz (1998) [8] (see also Natario (2000) [9]) in writing 
B k = rf^djAk, where rfi = (v/7) -1 e ■ The 4- vector B x is invariantly defined by 
B \ = -Ti^dnfa/Z 7 )^ = {0,B k ). Then B 2 = ±djA k V [j A k] where square brackets around 
indices denote the anti-symmetric part. Landau & Lifshitz [5] denote diAj — djAi by f\j so, 
translating their expression of Einstein's equations to our notation and writing £ = —V ln£, 

e a (r 1 v a e+ie 2 B 2 ) = i2oo, (15) 

- fW^ + -^ 2 {i kl B 2 - B k B l ) + P kl = R kl , (16) 

-^ 2 [Vx£-3£x£] fc = i? fc , (17) 

where P kl is the curvature tensor of the 7^ 3-space formed from 7^ just as R^ u is formed 
from in 4-space. Now R = g^R^ = {w^w u - ^^)R^ U = £~ 2 R o ~ lkiR k \ where 
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lnv = w^w v — g^u is a 4-dimensional covariant version of 7^. From (fl5j) and (|T6l) 

2rv 2 e - ^ 2 # 2 -p = r, (is) 

from (USD and (USD 



*T 00 = Poo - ^ = e {^ 2 B 2 + . (19) 

Now in the introduction we found that an important step in making the field energy a perfect 
square was the introduction of a conformal transformation which yielded a flat 3-space. We 
can not do that generally but we can transform so that our new 3-space has a vanishing 
scalar curvature. We write 7 fcZ = e~ 2 ®^ki and use the relationships between the curvature 
of two conformally related spaces given, e.g., in Stephani et al. (2003) [10], equation 3.85 
(contracted) 

e -2<Kp = p + 4V 2 $ _ 2 (v$) 2 . 
So the transformation that makes P zero obeys 

- 2V 2 $ + (V$) 2 = ^P, (20) 
which we may rewrite like a Schrodinger equation by setting ip = e~2* so that 

V 2 ^ - \p^ = . (21) 

o 

Equation (|20|) is the generalisation of equation ([9]) of the introduction and using (fl~9l takes 
the form 

- V 2 $ = ~ [r 2 ^ 00 - \i 2 B 2 - |V$| 2 j , (22) 

which shows how both gravomagnetic and gravitational field strengths subtract from the 
'source' of $. However previously that source was written in terms of kTq not kT 00 and 
since our expression (fl4l) for mechanical energy also involves kT$ we now rewrite our basic 
equations (1T5T) (TT8T) in terms of kTq , kTq and nT kl . To do this we use equations (fT9l) and 
f|T7|) and aim to get only second derivative terms on the left kTq = £~ 2 kT 00 + A k B,Q, so 

^ 2 B 2 + \p~ -(VxB-3SxB) = kT° . (23) 

Using (EQD for P and -\£ 2 A • V x B = \ V • {i 2 A X B) - \B ■ V X {£ 2 A) , we find on simplifying 
the right-hand side of (T5i 



V ■ xB- 2V$) = V • V = k(T ° + PG ), (24) 



where 



and from (TT3 



np G = -\eB 2 - | V$| 2 - \eA • (£ x S) ; (25) 



(V x Bf = -2kC 2 T* + 3(£ x B) k . (26) 



Notice that all quantities on the right of (l25j) are defined in terms of a Killing vector and 
(the square of the gradiant of) a conformal factor which, being the solution of the 'invari- 
ant' elliptic equation (120]) . depends only on the slicing, its geometry and on the boundary 
conditions. 

Equations ( |24l) and (26) may be thought of as a non-linear generalisation of Maxwell's 
equations and indeed we shall see presently that the terms of ( 1251) constitute the energy 
density of the gravitational field. 

The final equation of this trio comes from ffTB"]) and f|T8|) and reads 

r X (7 W V 2 £ - V fe V'0 = kT m - \^P kl - ]p H Pj ~ (]p kl B 2 ~ B k B l ^j . (27) 

Now in making our conformal transformation we recognise that our space metric will tend to 
the Schwarzschild form at infinity with m = GM being the Schwarzschild asymptotic mass 
parameter. Since Schwarzschild space is conformally 3-flat we may impose the boundary 
condition $ — > 0(l/r) on the solution $ of fl2"U|) . Then $ tends to the corresponding 
Schwarzschild value m/r found in equation ([8]) of the introduction as r —>■ oo, that is we 
identify the coefficient of 1/r as the total mass parameter. When we integrate equation (|24|) 
over all space so as to generate the mechanical energy Em from the first term on the right, 
we find on the left 

J V • xB- 2V$^) dV = j (]^ 2 A xB- 2V$^) ■ dS, 

where the dS integral is to be evaluated over the sphere at infinity still assuming no black 
hole is present. Now B is 0(l/r 3 ) there, the triad component of A will be 0(l/r 2 ) like the 
classical vector potential and £ 2 — > 1 so the integral of the first term on the right vanishes; 
however $ — > m/r, so the second term gives 8irm = kM. Using this result, the integration 
of ([21]) over all space yields 

Eg — M — Em — j p G dV , (28) 
10 



with pc given in (1251) . 

As KLB [6] pointed out the mechanical and gravitational energy densities defined above 
fail for systems with ergospheres because they involve = in their definition and 
£ is zero on the ergosphere. However such difficulties can be circumvented by using the 
mechanical energy as estimated by hypersurface orthogonal observers such as Bardeen's 
ZAMOs. We take such observers moving orthogonally to surfaces of constant time with 
4-velocities = C -1 (l) W k ). The appropriate form of metric is in terms of lapse and shift 
written as the final expression of (|T2j) . The metric components and their inverses were given 
in equations 2.37 and 2.38 of KLB [6] and we notice that Wk = 

9oo = ( 2 _ W 2 , g ol =W t = % t W k , g kl = -%i , 

g 00 = C -2 ^ g 0l = C 2 W 1 ^ g M = _~M + C 2 w k w l ; 



V~9 = (Vl , 7 = detjki- (29) 

The Einstein equations are 

kT 00 = G 00 = -C\P + K 2 - K kl K H ) , (30) 

kT* = G° k = -C'ViiKl - 6 l K K) , (31) 

-nT = R = 2(C 2 G 00 -f l Rki) , (32) 



kT u = Rki = Ph + V (fc £ - SA + 2(K$V [l) W m] - C% n W m K kl + V m (W m K k l) , (33) 

where Pki is the Ricci tensor of the spatial metric 7^, P = ^ kl Pki is the corresponding 
3-scalar curvature and Kki is the second fundamental form of the hypersurface t =constant. 
Thus 

Kki = C'V^) , K = K k = f l K u , and S k = -d k ln( . (34) 

The mechanical energy density on any hypersurface E of normal n M is for dust p(u.n) 2 
so for our observers moving orthogonally to the cut n fl = and for general T^ v the 
mechanical energy on an element of hypersurface yf—gdU^ is T^Wy^f—gdT,^ = C^T^yf^cPx. 
This is given directly by equation (J30l . Once again P is determined by using a conformal 
transformation $ such that the transformed 7^ has vanishing 3-curvature. Analogously to 
fl20l) this leads to 

-2V 2 4> + |V<f>| 2 = -P . (35) 
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Transporting second derivative terms to the left this yields in place of (!30|) 

K( ,2 T oa + 2 y2,| = |y^| 2 + _( R 2 _ K kl K kl ) . (36) 

Integrating and using the boundary condition that V$ — > (GM/r 2 )r at infinity we find for 
maximal slices (K = 0) 



E G = E M -M = - 

K 



|Vlf - l -K kl K k 



dV . (37) 



The mechanical energy density Em has real advantages over Em- Not only can it be mea- 
sured within the ergosphere but also the two '7' factors (u.n) are the same and correspond 
to what a hypersurface orthogonal observer would see. It also coincides with the geomet- 
rical expression of Wald alluded to earlier. However there are strong arguments against it 
also. The observers move relative to static observers and worse still they move relative to 
each other. As seen from infinity they circulate and for the axially symmetric case with 
azimuthal Killing vector 7/ M they rotate about the axis at angular velocity u = ^.77/ ( — T7.77) 
which depends on position. In this case the space-time cuts correspond to surfaces of con- 
stant Boyer-Lindquist time. Now in classical physics observers who rotate at angular velocity 
O see as their energy not the true energy E relative to observers at rest at infinity but rather 
the Jacobi constant E — Q.J, where J is the total angular momentum. For observers such 
as ZAMOs in differential rotation we expect E — J ujdV, where j is the angular momentum 
density about the symmetry axis. It is interesting that our expression for Em takes just 
such a form (writing dV = ^d^x): 

E M = J T^w^gd^ = J T^dV = J T?(C + cur)»)dV = T-J ujdV, (38) 

where T = J T^^/^jd^x is the part of the mechanical energy independent of the angular 
velocity of the observer and j = —T^nf is the angular momentum density. If, following this 
line of thought, we identify T with the mechanical energy rather than Em then we need to 
evaluate T® rather than T 00 . Einstein's equations in mixed form follow from (|30|) and (|3T|) 
by writing G® = (g 00 )" 1 ^ 00 — g 0k Gl); hence using ( 135]) and putting all second derivatives 
on the left 

KT*+2V 2 §-Vi\C 1 W k (K kl -^ kl K)} = \^\ 2 +^(K 2 -K kl K k i)-C 1 W k (K l k -6 l k K)ei. (39) 
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On integration the third term vanishes over the sphere at infinity and in the absence of black 
holes we have, setting K = for a maximal t =constant slice, 

-V = T- M = iy [|V$| 2 - \k u K u - C'WkK^S^dV, (40) 

which expresses the new gravitational potential energy V in terms of 'field' quantities. 
Advantages of T as a mechanical energy are: 

(i) like Em it can be evaluated for systems with ergospheres; 

(ii) like Em it involves mu • £ which is the energy for a dust particle, even for one that 
darts into and out from the ergosphere. 

(iii) It has removed the part of Em which is related to the circulation of the ZAMOs 
around the axis. 

However, even for dust it is not clear that it can be related to kinetic energy as seen by 
any chosen observers. 

III. MECHANICAL ANGULAR MOMENTUM DENSITIES 

For axially symmetrical systems there is no difficulty in defining total angular momentum. 
It is — J T£r] v ^—g dE^ and this translates into — J T^rj u (dV, or, for stationary systems 
without ergospheres — J Tjfrf^dV . However, whereas the split between mechanical energy 
and gravitational energy was clear, the split between mechanical angular momentum and 
field angular momentum caused us difficulty Luckily considerable enlightenment comes 
from first studying the electrodynamic analogue in flat space, so to this we now turn. We 
consider a charged rotating fluid held together by some cohesive force such as surface tension; 
a charged oil drop might be a good example. We wish to split the total angular momentum 
into a part due to the mechanics of the fluid itself and an electromagnetic part. We start by 
considering a single particle of mass m and charge q. Its 3-velocity will be v or in cylindrical 
polar coordinates (R,R(p,z), but that is written in local triad language; a relativist would 
use the 3-metric da 2 = dR 2 + ^^dcf 2 + dz 2 and write the 3-velocity as v k = (R, 0, z) or Vk = 
(R, 7000, i), where in this case 7^ is just R 2 . As recorded by static observers the mechanical 
angular momentum of a particle is just rwy^dcf) / dr) , where dr = dtyl — v 2 ; however this is 
not the total angular momentum associated with the particle because it is charged and moves 
in an electromagnetic field. Its Lagrangian will be — my/1 — v 2 + qv. A — qA where A is 



13 



the electromagnetic vector potential and Aq is the electrostatic potential. Thus dL/dcf) = p^ 
will have two pieces: rwy^fi/y/l — v 2 + qA^. The first is the mechanical angular momentum 
m l<f>^>(d(f) I dr) , considered above; the second is a piece of electromagnetic angular momentum 
that is to some extent associated with the motion of the particle and its charge, but it also 
depends on all the other charges that act together to give rise to the vector potential A. 
Notice this piece of momentum is not generally gauge invariant because we have not yet said 
anything about fixing the electromagnetic gauge. It is actually incorrect to believe that the 
sum over all particles of all their p^ gives the total angular momentum. It does not! The 
trouble arises just because the A^ is generated by other particles in the same assemblage 
- we encounter a similar problem when adding individual particle energies in which the 
electrical potential is included. Even classically the straight sum counts the electrostatic 
energy twice. What then is the correct procedure in the electrodynamic case? The answer 
is to separate the mechanical angular momentum which is additive and gauge invariant. We 
sum this over all the particles. Then quite separately we work out the total electromagnetic 
angular momentum from the gauge invariant expression J em = j[r x (Ex H)].zci 3 a; = 
j M^rj^d 3 x, where M/f is the Maxwell stress-energy- momentum of the electromagnetic field 
and rf is the angular Killing vector, so — 1. When we sum Jm + J em we find that we do 
indeed get the total angular momentum. This is obvious since the total T^ v is the sum of 
the mechanical and the electromagnetic energy momentum tensors. This final more obvious 
option is not open to us when we deal with the gravitational field's angular momentum, as 
there is no known stress energy tensor for it; however, as we now show, the earlier argument 
for splitting the mechanical angular momentum from the field part can be followed in the 
gravitational case to which we now return. 

We shall again consider a single particle but now it will be uncharged and in a stationary 
metric as in ffl2|) . i.e., 

ds 2 = i 2 (dt - A k dx k ) 2 - -y k idx k dx l , (41) 

but we shall take this metric to be axially symmetric and, to start with, we consider the 
still simpler case in which there is a <fi — > —<p, t —>■ —t symmetry. In the latter case we write 
Akdx k = A^dcj) + A^dx 1 and 

r )tidx k dx l = 700<i0 2 + jKLdx K dx L , (42) 

where K and L run from 2 to 3 whereas k, I run from 1 to 3 and Greek suffices run from 

14 



to 3. 



IV. MECHANICAL ANGULAR MOMENTUM FROM STATIC OBSERVERS 

To keep contact with directly observable quantities we introduce our set of static observers 
(see Section II). The speed of the fluid, v, past the observers is given as (1 — t> 2 ) -1 / 2 = 
w^u^ = (w.u) and the components of yp=f are given by u — (u.w)w. Hence the mechanical 
momentum of the particle, unconnected with its field momentum, will be m[u— (u.w)w] 
and the corresponding mechanical angular momentum is m(u — (u.w)w).?7, r\ is the axial 
Killing vector. Since only the component of r\ transverse to £ M is involved in this product, 
this angular momentum can also be re-expressed as m u^^^rf , where 7 Mi/ = w^w u — g^ v . 
If we consider a stress-energy tensor made up of dust, then T^ v = pu^u u and pu^ is the 
rest mass flux vector; multiply by u u, y ucr r] cr , the mechanical angular momentum per unit rest 
mass, and pu fJ, u u, y ucr r] cr = T fJ-u, y ucr r] cr is the flux of mechanical angular momentum. We shall 
now find this by another method directly analogous to our electrodynamic calculation. 

The Lagrangian for a free particle moving in the metric (|4T|) is L = — m[£ 2 (l — A^<p ~ 
Akx k ) 2 — l4>4>4> 2 — 1kl% R x L } l l 2 = —m ds/dt. Its total angular momentum will be p^ = 
dL/dcj), so 

P<j> = m[£ 2 (l - - A K x K )A^ + ^ H 4>]/ (ds/ dt) , 
= m£ 2 (u° - A+u* - A K u K )A$ + rrvy^u* . 

Now u» = (dt/ds, dx k /ds), so u = £ 2 (u° — A^vft — Aku k ) and {l-v 2 )' 1 / 2 = w^u^ = f _1 u , 
so reversing the order of the terms we get p^ = m/y^u^ + m£ (1 — v 2 )~ 1//2 A/>- This expression 
is in precise correspondence with the electromagnetic case with the mechanical angular 
momentum in the first term and the gauge dependent field term coming afterwards. When 
the system does not have — > —(f), t — > — t as a symmetry, there are dcf)dx K terms in 
r Ykidx k dx l , and the form of p^/m becomes 

P4>/m = ^ k u k + f (1 - v 2 )- 1/2 A^ . 

The mechanical angular momentum per unit mass is the first term and to get the angular 
momentum flux vector we multiply this by the rest mass flux pu 11 , so the mechanical angular 
momentum flux remains T^j^ = T^j^ = T^ v ^ va if . 
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Although for easy explanation we have adopted the simplest dust case to explain our 
points, nevertheless our final formulae still hold whatever the constitution of T^ v . 

The total gravitational field angular momentum reckoned by static observers is the dif- 
ference between the total angular momentum and the total mechanical angular momentum, 
so for a spacelike surface S 

Jgs = J-Jms = j -r^^+v)^^^ = J -T^WvWeify/^gdZn = - J T^A^dV 

Now we have established this as the angular momentum they attribute to the gravitational 
field our next aim is to re-express it in terms of the squares of field quantities by using 
Einstein's equations for T^ u and performing integrations by parts. Multiplying (124)) by ^A^, 
integrating over all space we find on smuggling the ^A^ into the divergence and then paying 
the duty 

J ZAP -dS-J V(£A ) • VdV = k(J T^A^dV + j p G ^dV 

Now A<p is 0(l/r) and T> is 0(l/r 2 ), so the first term vanishes when integrated over the 
sphere at infinity while the first term on the right gives — Jgs- Hence, 



J 



GS 



p G £A + -V(^) • V 



dV, 



where V is given in (EHj) and pc in (125|) . When A has only a <fi component, VkA^ = —r\^ m B m , 
so then V ■ V0 = — (T> x B)<j, and the final term becomes -[—(T> x B)^ — T>- EA<j\. 



V. MECHANICAL ANGULAR MOMENTUM WITH RESPECT TO ZAMO'S 

The ZAMOs have a different concept of what constitutes space than the static observers 
and that concept is superior in that it extends within ergospheres. It is natural to replace 
the w u of the static observers in equation (TH1) by the w u = n v of the ZAMOs. Indeed if we 
calculate in the appropriate ZAMOs coordinates, we find 

P4> = ^- = m% k (u k W k u°), 
ocp 

and it is just the u k term that the ZAMOs consider to be the mechanical angular momentum 
with respect to them. Their estimate of mechanical angular momentum is therefore (with 
(p unsummed) 

Jm = J T^~ lva r?^r g d^ . 
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But J = — J T^ h 'g U aV° V—QdT,^ and rf has only a component and since 7^ = —g^v+w^Wu, 
it follows that Jm — J = f T^ v w v w a rf y/—g dE^ = because w a rf = 0. 

Thus the total angular momentum is the same as the material angular momentum as 
seen in the local rest frames of the ZAMOs. 

The mechanical angular momentum defined via static observers cannot be assessed for 
systems with ergospheres. The mechanical angular momentum defined via ZAMOs takes 
no account of their circulation around the axis. Is there a splitting into mechanical and 
field angular momentum that takes account of the ZAMOs motion but does not break down 
within ergospheres? Something like the split we found for energy by writing = £ M + cur/ 4 
is needed. 

For a dust fluid the angular momentum is J = — J pu^u v rf \J^g dS M = — J pu°u u ' g v $dV ' . 
Now u u g V! p = u°g 0( j, + u k g^ = u°W<p — u k ^k<p- The first term on the right is due to the 
circulation of the ZAMOs, the last is most readily seen as an angular momentum per unit 
mass of the fluid when u k has only a 0-component. Thus, following Bardeen (1970) [7], J 
can be split as follows 



For those who want this split expressed covariantly the first 'field' term may be rewritten 



To evaluate the field angular momentum we rewrite (36) in a form suitable for construct- 
ing the required integral (setting K 2 = for a maximal slice): 



where jf is defined in equation fH3l) . ( and can be expressed in terms of the Killing 
vectors £ and the angular one r\ via — £ • T] and ( 2 = £ 2 — (W^/rj 2 ). 



J 






Jf 
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VI. CONCLUSIONS 



We have elucidated the physics of three different expressions for the energy and angular 
momentum densities of gravitational fields. While the authors differ as to the relative 
importance of these expressions DL-B favours the third. 

None of these densities make contributions to T^ v although they do act as (negative) 
sources for the conformal field $ which gives the scalar curvature of 3-space. 

The fact that they do not contribute to the source T^ u was one of MTWs strongest 
arguments for the non-existence of a stress-energy tensor for gravity. That argument still 
stands. We have conducted our investigation within General Relativity. A different theory 
might have the field energy contributing to the source T^ u but that is not General Relativity. 
Energy densities are only really useful when the system changes, whereas this paper is 
confined to stationary systems. Nevertheless, we believe these methods may be generalised 
to systems that change so slowly that no gravitational waves are emitted. 
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